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ABSTRACT 

A systematic method is given for obtaining consistent approximations to the 
Schwinger-Dyson(SD) and Bethe-Salpeter(BS) equations which maintain the ex- 
ternal gauge invariance. We show that for any order of approximation to the SD 
equation there is a corresponding approximation to the BS equations such that the 
solutions to those equations satisfy the Ward-Takahashi identities of the external 
gauge symmetry. This formulation also clarifies the way how we can calculate the 
Green functions of current operators in a consistent manner with the gauge invari- 
ance and the axial anomaly. We show which type of diagrams for the tx^ 77 
amplitude using the pion BS amplitude give result consistent with the low-energy 
theorem. An interesting phenomenon is observed in the ladder approximation that 
the low energy theorem is saturated by the zeroth order terms in the external 
momenta of the pseudoscalar BS amplitude and the vector vertex functions. 
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1. Introduction 



If we consider a strongly interacting fermion system, we have to deal with 
various boundstates and their interactions among themselves or to some external 
gauge fields. For instance, in QCD, the boundstates are hadrons, which are in fact 
the only observable states, and the external gauge fields are the photon and the 
W and Z bosons. In order to treat such boundstate problems, we are obliged to 
adopt some non-perturbative approximation scheme, and then it becomes a non- 
trivial issue whether and how to keep the external gauge invariance. 

In this paper we discuss this problem for the case using the Schwinger-Dyson 
(SD) and Bethe-Salpeter (BS) equations as a non-perturbative approximation 
method [1]. We use, in a certain approximate manner, the SD equation for cal- 
culating fermion propagator and the inhomogeneous BS equations for the vertices 
of fermion with the external gauge bosons. Then the adopted approximation has 
to keep the external gauge invariance in order to give results consistent with, for 
instance, various low-energy theorems. This problem is actually not so trivial and, 
in fact, it will turn out that the approximations for the SD and BS equations can- 
not be independent of one another but have to satisfy a mutual consistency. This 
is because the fermion propagator and the vertex functions are connected to one 
another by the Ward-Takahashi (WT) identities resulting from the external gauge 
invariance. We shall present a general answer to this problem and find that for any 
order of approximation to the SD equation there is a corresponding approximation 
to the BS equations which satisfies the external gauge invariance. 

A partial answer to this problem has actually been known to Maskawa and 
Nakajima [2, 3] for a long time: the axial- vector vertex function satisfies the axial- 
vector WT identity if it is determined by the ladder BS equation in which the 
fermion propagator determined by the ladder SD equation is used. They proved 
this result by a direct diagrammatic evaluation of the divergence of the axial- vector 
vertex in that approximation. Our result in this paper is to give a generalization 
to arbitrary order of approximations. The proof of the gauge invariance is given 
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in a very systematic manner. 

Gauge invariance gives constraints not only on those fermion propagator and 
vertex functions but also the Green functions of the vector or axial- vector current 
operators coupled to the external gauge fields. We also discuss such Green functions 
of current operators in the same framework. We consider, for instance, a 3-point 
Green function of one axial-vector current and two vector currents, which gives 
a well known anomaly and is related the 7r° 77 amplitude by the low-energy 
theorem [4]. In the ladder approximation, as we show, the Green function simply 
obtained by a triangle diagram in which the vertices and fermion propagators are 
those in that approximation, satisfies the vector gauge invariance and reproduces 
the correct axial anomaly. This simple situation is, however, no longer true in 
any approximation beyond the ladder. We shall show that we have to include 
corrections intrinsic to the Green function which can be attributed neither to the 
fermion propagator nor to the vertex functions. This is also the case, even in ladder 
approximation, for general n-point Green function of current operators with n > A. 
For instance, consider the 4-point Green function of vector currents corresponding 
to photon scattering 77 — 77. The simple box diagram consisting of the vertices 
and fermion propagators calculated in the ladder approximation does not satisfy 
gauge invariance, despite that the vertex function and the propagator themselves 
are mutually consistent with WT identity in the ladder approximation. We shall 
clarify in a general manner which types of diagrams should be included in order to 
obtain gauge- invariant Green functions (with correct anomaly in case it is present) . 

We discuss, in particular, the calculation of 7r° 77 decay amplitude by using 
pseudoscalar BS amplitude in detail, since it gives a typical example in which 
all these problems of gauge invariance and anomaly become relevant. We show 
that the resultant amplitude satisfies the low-energy theorem if the boundstate 
BS amplitude is the one determined by the BS equation in the 'same' order of 
approximation as those for the fermion propagator and the vector vertices, and if 
we include a particular set of diagrams which depend on the approximation. We 
also point out an interesting phenomenon in the ladder approximation that the low 



3 



energy theorem is saturated by the zeroth order terms in the external momenta of 
the vector vertices and the pseudoscalar BS amphtude. 

This paper is organized as follows. In section 2, we introduce a gauge invariant 
effective action for the fermion propagator on arbitrary background of external 
gauge fields. This provides us with a very systematic way to obtain mutually 
consistent SD and BS equations. It is proved that the solutions to these equations 
satisfy the WT identities in section 3. In section 4, we construct 3-point Green 
function of one axial-vector and two vector current operators, and show that it 
not only satisfies the vector gauge invariance but also reproduces correct axial 
anomaly. Section 5 is devoted to consideration of the pseudoscalar boundstate 
and the low-energy theorem for 7r° — > 77. We show which types of diagrams give 
amplitude consistent with the low-energy theorem. In section 6, we discuss the 
above mentioned phenomenon of zeroth order saturation of the low-energy theorem 
in the ladder approximation. An appendix is added to show how the discussion of 
the anomaly in section 4 goes when we adopt dimensional regularization instead 
of the Pauli-Villars-Gupta's one adopted in the text. 

2. Gauge Invariant Effective Action 

We consider an interacting fermion system like QCD in which the chiral sym- 
metry is spontaneously broken dynamically. As usual we call the 'gauge' interaction 
responsible for the formation of boundstates and for the spontaneous chiral sym- 
metry breaking "color gauge interaction" and the gauge boson "gluon" , although 
the present formulation also applies to more general systems than QCD.* In such 
a system we want to calculate, in an approximate but non-perturbative manner, the 
fermion propagator as well as the (color-singlet) vector and axial- vector vertices 
to which external gauge fields couple. A method is to use the Schwinger-Dyson 

★ The 'gauge' boson need not be a true gauge boson; so, for example, it may have non- 
zero masses and may be axial- vector. Then the Nambu-Jona-Lasinio [5] like models can be 
discussed by considering a limit the vector and axial- vector 'gauge' boson masses go to very 
large. 
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(SD) equation for the fermion propagator and the inhomogeneous Bethe-Salpeter 
(BS) equation for the vertices. In this announced in the Introduction, the 

approximations adopted for the SD and BS equations cannot be independent one 
another but have to satisfy a mutual consistency in order to meet the the external 
gauge invariance requirement. 

A systematic way to obtain those SD and BS equations satisfying the mutual 
consistency, is provided if we consider the original system put in a general external 
gauge field background. Let us denote the external background gauge fields as 

= A^^y^ with fiavor matrix A'* normahzed as tr(A'*A^) = {l/2)6ab: ^-nd assume 
a vector-like couphng to the fermion: >Cint = 'tp^nA^^ip. We assume this just for 
notational simplicity and axial- vector case can be obtained simply by replacement 
7^ 7^75- An important assumption is that the flavor degrees of freedom to which 
the external gauge fields couple arc orthogonal to the color degrees of freedom to 
which the 'gluons' (internal 'gauge' bosons) couple. 

Effective action for fermion propagator Sp in the presence of external gauge 
field background is given by [6, 7] 

r[SF,A] = iTrLn,5F - Tr (i^-Sp) + r^/C2Pi[-SF], (2.1) 

where we note that the external gauge field A^ is present only at the covariant 
derivative — — iAj^. Here /C2P1 stands for the two particle irreducible (w.r.t. 
fermion-line) diagram contributions: for example in QCD, we can expand the /C2P1 
into power series of the gauge coupling as = 5's/47r, 

/C2PI = /Cgi + /C^p\ + • • • (2.2) 

and /Cgpj and /Cgpj are diagrammatically given by Fig. 1. 
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Fig. 1. Two particle irreducible (w.r.t. fermion-line) diagrams 
contributing to IC^Pl ^2PI ™ QCD. The double wavy 
line represents the gluon propagator D^, and the solid line 
represents the fermion propagator Sp. 

More exphcitly the first term /Cgpj is given by^ 

2 p 

4pi = -f y '^'^'^'^ l/)nMr"^F(t/, x)z7.r«) ^'^^(x - t/) , (2.3) 

where (a = 1,- ■ - .iV^ _ 

arc color matrices in the quark representation and 
D^i, is tree level gluon propagator given by 

(2.4) 

where we have included an ultraviolet cutoff A for definitencss. For later con- 
venience, we refer to the first diagram contribution to the K-^Pl ^^S- ^0^) ^ 



f If we use the running coupling constant, the couphng constant gs depends on the gluon 
momentum and hence we should understand that the running coupling factors in such 
a case are placed at the loop integrand. All the discussions in this paper apply equally 
both to the fixed and running coupling cases, provided that in the latter case the gluon 
momentum is used as the argument of the running coupling function [3]. 
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j^(2a) 

X tr (^5'F(xi,|/i)i7^T"5F(i/i,X2)i7pT''5F(x2,y2)«7i'7'"5'F(i/2,a;i)i7aT^' 

(2.5) 

By our assumption that the flavor is a freedom orthogonal to the color, we note 
that the flavor matrices A'* commute with the color matrices T". Then we have 
the following lemma: 

Lemma: For any approximation for /C2P1 by an arbitrary subset of diagrams con- 
tributing to /C2PI7 the effective action Eq.{2.1) is (external) gauge invariant: 

r[5F,A]=r[5F,^^] , (2.6) 

where the gauge transformation with U{x) — exp {i9°'{x)\'^^ is given explicitly by 

SF{x,y) ^ S^{x,y) = U{x)SF{x,y)U-\y) . 

Proof) It is convenient to introduce the following functional notation: 

iU)^y ^ U{x)8\x - y) , 

{p),y^ip.-iM^))^\^-y). (2.8) 

In this notation we can rewrite the gauge transformation (2.7) simply in the form 

p upu-'^ , 

(2.9) 

Sf USfU-^ . 

Then we can see the gauge invariance of the first two terms in F as follows: 

TrLn-Sp ^ TrLn {USfU~'^) = Tr [C/ (Ln^p) U''^] = TrLn-Sp , (2.10) 
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Tr (^^f) ^ Tr {UpU-^USpU"^) = Tr {pSw) • (2.11) 

The gauge invariance of /C2Pi[>S'f] is also shown similarly: for instance, for the 
lowest order diagram we have 

4pi^ -f / d'^d^y {U{x)SF{x,y)U-\y)i^^T^U{y)SF{y,x)U-\x)i^,T^) 
xD'^''{x-y)=}C^^,. 

(2.12) 

The point here is that the fermion propagators appear successively in the trace and 
the vertices there contain only color and 7 matrices which commute with the flavor 

gauge transformation matrices U{x). This property clearly holds for any diagrams 
contributing to /C2P1 and so the gauge invariance follows. This finishes the proof. 

Because of the lemma, the Schwinger-Dyson (SD) equation derived from the 
action F is automatically (external) gauge covariant (or gauge invariant as a set of 
equations) . 

The SD equation is given by dT/dSp = 0, which reads 

^S^' = ^P-^-'^-^. (2.13) 

If we take only the lowest order term in /C2P1, ^2PI' then this SD equation reduces 
to (see Fig. 2) 

iSp^ = i0 +^ + r^K *Sf , (2.14) 

with K * Sf defined by 

K*SF = gUt7^.T^)SF{y,x){t^,T'')D^''{x-y) . (2.15) 



Eq.(2.13) (or (2.14)) is the SD equation determining a solution Sf — Sf[A] for the 
fermion propagator, on an arbitrary external background gauge field A^. 
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Fig. 2. Schwinger Dyson equation derived from the effective action 
r using /C2P1 = /C2Pr 

Tfie solution S-p[A] is expanded into a power series in the external gauge field 

A, 



Sf[A] =Sf + lA^Gf + '^A'l^AlGf'^'' + '^AlAlApf^^^'^^P + • • • , (2.16) 

where a, h and c denote the flavor indices. Here and henceforth the space-time coor- 
dinates and the integrations are suppressed, i.e., A'^G'^^ = J d'^zA1j^{z)G'^'^{x, y; z), 
etc.. The function G^l^g' "'"'"'^"(x, y; 2^1, • • • , Zn) defines a fermion 2-point function 
with n vector vertices inserted (see Fig. 3): 



G,^(x,y;z)=- ^^^^^^ 



^{0\Tf^{z)i^{x)i;{y)\0) , 

A=0 



G^'^'^^^fa; vzw)^ 1 SSFix,y;A) 



= {0\Tf^{z)j'''{w)ij{x)ij{y)\0), .... 

A=0 

(2.17) 

This is because 5/5A1^ yields an insertion of the vector current operator j'*^ = 
■07^A'^'0 to which the external gauge boson couples. Hereafter in this section, 
we suppress the flavor indices to denote (^^i^^^'"'-'*"^" simply as ^^+2^^", and write 
only in place of 7^A'* as vertex factors in the figures, accordingly. 

Therefore the functional differentiation w.r.t. (and then setting A = 0) 
of the SD eq.(2.14) successively generates the Bethe-Salpeter (BS) equations for 
the functions and they are automatically (external) gauge covariant. It is 
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a) 



b) 



Fig. 3. Graphical representations of a) Gg and b) G'^^ defined in 
(2.17), where wavy hne represents the external gauge field. 

convenient to define the following vertex function by amputating the fermion legs: 



^ n+2 — "-"f '~^F ■ 



(2.18) 



To show what is going on as explicitly as possible, from here on in this section, we 
confine ourselves to the simplest case using the lowest order kernel (/C2PI = ^2Pl)- 
The extension to more complicated case will be trivial. First differentiation S/SAfj, 
of (2.14) gives (see Fig. 4) 



(2.19) 



where K*r'^ = X * (^'pFg^'p) = K*G'^ is defined in the same way as in Eq.(2.15). 





+ 




Fig. 4. BS equation for T^. 
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Fig. 5. BS equation for 



Second differentiation of (2.14) gives (see Fig. 5) 

rr - Tg^pr^ - r^^pF^ = k * Ff . (2.20) 

Third differentiation of (2.14) gives 

Ff ^ - (Ff >5fF^ + F^-SpFf + ill, u, p permutations)) 
+ (F^,SfF^5pF^' + (/X, I/, p permutations)) (2.21) 
= k * Ff ^ . 

We can see by the help of Eq.(2.20) that the last terms of the LHS in Eq.(2.21) 
play the role of cancelling the double counting of some diagrams contained in the 
second terms. Thus this equation reduces to 

Ff ^ - ({k * Ff )5'fF^ + F^5f(X * F^'') + (/X, I/, p permutations)) 
- (F^-SpF^-SpFf + (//, V, p permutations)) (2.22) 
= i^*Ff^ . 

Eqs.(2.19), (2.20) and (2.22) are the BS equations determining the vertices 
Fg , F^^ and Fg'^'^, respectively, in the ladder approximation. It is in fact easy 
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Fig. 6. Definition of L and L' . 

to find formal solutions to these equations. Define the four-point fermion Green 
function L in the ladder approximation as given in Fig. 6. 

Then, inspection of the Eqs.(2.19), (2.20) and (2.22) shows that the formal solutions 
for Fg , T'^ and Fg'^^ are given by Figs. 7, 8 and 9, respectively. 




' — V 

Fig. 7. Formal solution of Fg using X^2Pl — 
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Fig. 8. Formal solution of F^'^ using /C2P1 = ^2Pr 
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Fig. 9. Formal solution of Fg^^ using /C2P1 = ^2Pr 
3. Ward-Takahashi Identity 

We now show explicitly that the (external) gauge invariance of the effective 
action (or the covariance of our SD and BS equations) implies that the vertex 
functions determined by those BS equations satisfy Ward-Takahashi identities. 



Gauge invariance of the effective action F[S'f, A] implies 



r[SF,A]^r[s^,A^] 
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(3.1) 



with = USpU ^, which leads to the gauge covariance of the SD equation: 



SS' 



(3.2) 



Therefore, if Sp — Sf[A] is a solution to the SD equation on the background A, 
i.e., 

ST[Sf,A] 



6Sf 



Sf=Sf[A] 



0, 



(3.3) 



then, its gauge-transformed one = USp[A]U ^ gives a solution to the SD 
equation on the gauge-transformed background A^: 



ST[S^,AU] 



. 



(3.4) 



S^=USf[A]U-'' 



Namely, we have shown 



Sf[A^] = USf[A]U-'^ . 



(3.5) 



If there are several solutions Sp[A] {i — 1,2, • • •) for a single background A, we 
should have relation ^'pfA'^] = USp[A]U~'^ for each i, because of the continuity for 

Substituting the expansion (2.16) into both sides of (3.5), we have 



LHS = Sf + tAp^, + '-A^A^iG^^ + • • • , 

RHS = USfU^^ + iA^^UGlU-^ + -A^AJJC^^IJ-^ + • • • . 



(3.6) 



In particular, for an infinitesimal gauge transformation U — \^iQ {Q — 9°-\°-), for 
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which is given by 



(3.7) 



with structure constant fate of the flavor group, we flnd 



LHS = {Sf + zd^e'^G^'^) + [{Sab + facbOl Gf + idj'Gl^''''\ + 
RHS = ,5f + i {OSf - SyB) + iAl [Gl^ + i {OGf - GfO)] + • • • , 



(3.8) 



are equal with each other. Equating each power term in on both sides, we flnd 



(x, y- z) = i5\z - x)\^Sy{x - y) - i5\z - y)SF{x - y)\\ (3.9) 

-id^Gl^'^^ix, y; z, w) = -fabc^^iw - z)Gl^{x, y; z) 



. . ., and so on. These arc just the Ward-Takahashi identities required by the 
external gauge invariance. For instance, Eq.(3.9) is nothing but the WT identity: 



dl (0| T3-^{z)i,{x)i,{y) |0) = -S\z-x)\- (0| Ti,{x)i,{y) n^S\z-y) (0| T^(a;)^(|/) |0) A'^ .| 



Thus this proves that the fermion propagator Sf and the vertices F^^g'^" deter- 
mined by our SD and BS equations satisfy the Ward-Takahashi identities giving 
relations among them; namely, our approximations for the SD and BS equations 
are mutually consistent and gauge invariant. 

Note that we have proven that this external gauge invariance holds not only for 
the simplest ladder case (i.e., with /C2P1), but also for any order of approximations.* 
For example, if we take /C2P1 = ^2PI + ^2PI^' then the SD equation for Sf and the 
BS equation for Fg'^ are changed into the forms as shown in Figs. 10 and 11. These 
equations are much more complicated than the simple ladder ones, nevertheless 
they satisfy the gauge invariance. Important is the mutual consistency of the 
approximations between the SD equation and BS equations. 

■k The gauge invariance for the 3-point vertex F"^ in the simplest ladder case has been 
known for a long time to Maskawa and Nakajima [2]. A refinement of the proof and the 
genelarization to the running coupling case was given by Kugo and Mitchard [3] . 



-|- i6'^{w — X 



■)A^Gf (a:, y; z) - tS\w - y)G''/{x, y; z^X\0) 
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4. Anomaly 

We now turn to the problem of triangle anomaly. For definiteness, we work 

in SU{3)c QCD with two flavors; wc consider SU{2)i^ x 5*^7(2)^ chiral symmetry 
limit in which the u and d quarks have zero masses. 

Consider fermion 2-point function with a vector and an axial-vector current 
inserted, 0^^'^°^, where Q/j, denote the vector current to which photon couples, 

and ira denotes the axial-vector current to which pion tt^ couples: 
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Then the WT identity Eq.(3.10), which was proven for the vector case in the 
previous section and can straightforwardly be extended to the axial-vector case, 
implies the following vector and axial-vector WT identity: 



= i5\y -w)Q GTiw, z- x) - i5\y - z) GTiw, z;x)Q, (4.1) 
-td-,G'^^'^^{w,z;y,x) 

^ iS\x - w) |75Gf z; y) + i5\x - z)G^^{w, z- y) ^75 , (4.2) 

where the fabc term does not appear because of [Q, = 0. Therefore, if we close 
the open fermion legs of Q^i^''^'^ by inserting another vector vertex factor Qj'^, we 
obtain formally both vector and axial-vector gauge invariant 3-point function: 

T^^'^ix, y, z) = t^e^ tr (g'J^'""(^, z; y, x)QY) 

= i'e' tr ((0| Tf(y)ji(x)^(z)^(z) |0) Q^'') (4.3) 
= {0\Tji{x)fiy)fiz)\0) . 

Indeed the vector and axial- vector gauge invariance of T°^^^'^{x, y,z) is immediately 
seen from Eqs.(4.1) and (4.2): 

-idyT^^'''{x,y,z) 

= e^5\y - z) [tr {{QGf{z, z- x) - GI^{z, z- x)Q] Q^^)] 

= , (4.4) 

-idlT^>'-{x,y,z) 

= eH\x - z) [tr ({y75G^(^, z- y) + G^(^, z- yf-^^^] Q^^) 

= {formally) . (4.5) 



We should note two things here: 
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i) Despite its asymmetric looking in the above definition of T"'^'^, the two vector 
vertices Qfi and Qu are in fact on the same footing. Indeed, if the lowest order 
/Cjpi is used, for instance, recall that the formal solution for r4 was given in 
the form of Fig. 8. Taking it into account also that r3 was given by Fig. 7, we 
see that the amphtude T"^*^ is written in the following manifestly symmetric 
form with respect to the two vector vertices and stands for "triangle" diagram 
as shown in Fig. 12: 

rpa^„. = i^e'^Tr (^Fr'^^pr^^pr'^) + cross term) . (4.6) 

Here and henceforth, the vertices and F^'^ are simply denoted as Fg and 
F'^, respectively, for notational simplicity. 




Inspection shows that the same is also true for more general /C2P1; for exam- 
ple, when we take /C2PI = ^2PI ^2Pl''' formal solution for T"^^'^ take 
the form as given in Fig. 13, again showing manifest symmetry between the 
two vector vertices. 

ii) The above 'proof of the vector and axial-vector gauge invariance is just for- 
mal. The three point function T"^^ given in Eq.(4.3) is not well-defined as it 
stands: the Green function G4{w, z; y, x) itself is well-defined but (^4(2;, z; y, x)| 
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at coincident point w = appearing after closing the fermion legs, is a di- 
vergent quantity and no longer well-defined. 

The latter point is of course the well-known problem of anomaly. In order 
to define the three point function T°'^^^ properly, we need a regularization. Let 
us here adopt the Pauli-Villars- Gupta (PVG) regularization. (A brief discussion 
for the dimensional regularization is given in Appendix.) Since this regularization 
maintains the vector gauge-invariance, we have only to examine the axial-vector 
WT identity: QaT'^^^ —7. Hereafter we work with taking only the lowest order 
/C2P1 = ^2Pl' simplicity. Working in momentum space, we denote incoming 
momentum into the axial- vector vertex as qa, and out going momenta from the 
vector vertices and F^ as Pfj_ and kiy, respectively. In PVG regularization, the 
regularized 3-point function TZ'^\ is defined by 
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rtr 



(4.7) 



+ (cross term) , 

where we have omitted the obvious momentum arguments of the vertices and S-p 
is the propagator of the PVG regulator field with mass M: 



/-M 



(4.8) 



To see that this is really regularized, we actually need the high momentum behavior 
of the solutions Sp, F^, and . Since the gluon propagator is well regularized 
as stated in Eq.(2.4), the loop integrals in the SD and BS equations are sufficiently 
convergent. Prom those equations one can convince oneself by a simple power 
counting argument that the solutions approach the free ones as i — > oo: 



^f(^) 



ra 
5 



/+o(ri) , 



(4.9) 



T3 



7^75 + 0(^2) 



Therefore Eq.(4.7) gives a well regularized quantity: the degree of superficial diver- 
gence of the loop integral in this expression (4.7) is now only logarithmic (actually, 
the integral is convergent) and therefore the integral is independent of the shift 
of the loop momentum i [8] , contrary to the unregularized case. For definiteness, 
however, we choose the momentum assignment for the /i u cross term as shown 
in Fig. 12. Then, we can rewrite Eq.(4.7) into the form 



(Gf''^«(£,£-A;;p,?) 



- 5f(^)^7"755f(^ - q)QrSF{t - k))Qf 



(4.10) 
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with functional trace Tr implying also the momentum integration J d^£/(2TT)^. 
Using the momentum space version of Eq.(4.2) 

qa G?'^'""!^, i-k;p,q)= .^75 (£ -q,i-k)+ G^^{i, i + p)t^^, 

= i^^5SF{i - q)ri'SF{e -k) + SF{e)ri'SF{e+p)i^^5 

(4.11) 

and an algebraic identity 

^75 = S^\£ + p)i^5 + il^S^\^ -k) + 2M75 , (4.12) 

we find (see Fig. 14) 
^aT^Z^iP^k) 

= i^e^Tr 



.T3 

'2 



+ 



755f(^ - q)r''SF{£ -k)- z^75.5f(£ - q)Ql^SF{£ - k))Qf 
(^SF{i)r''SF{i+p)i^^5-SMQl^SF{i + p)t^75)Q7' 



- t'e^ lim 2MTrf 5755f(^ - k)-!" SF{t)l^SF{t + p) + (cross term)) . 

(4.13) 

The first trace term in the RHS vanishes if we can shift the integration variable 
i into f + g for the two terms in the first line since, then, they take the same form 
as the two terms in the second line and cancel them exactly by {2^75, (^7'^} = 0. 
This shift of the integration variable separately for the first and the second term 
yields nonvanishing surface terms [8], which are, however, the same for both terms 
and cancel each other thanks to the asymptotic behavior (4.9) of our solutions. 
The remaining last term reproduces the well-known anomaly: 

?«^?egV^' ^) = ^'i^^' ^^3QQ) e^-P'^ppk, = [Anomaly] . (4.14) 

We thus have shown that our 3-point function calculated by using the non- 
perturbative fermion propagator Sf and vertices and F^, not only satisfy the 
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Fig. 14. Graphical representation of axial-vector Ward-Takahashi 
identity, qaT^^^^{p, k), where dashed lines represent the 
insertion of i^Ts the place where the momentum q 
flows in, and the double solid lines represent the propa- 
gator of PVG regulator field, Sp. 

conservation for the vector channels but does in fact reproduce correct anomaly 
for the axial-vector channel, the latter being in accord with Adler-Bardeen theo- 
rem [9] . Note again that this nontrivial property is achieved only when the fermion 
propagator and the vertices are determined by such SD and BS equations with ap- 
proximations mutually consistent with each other. 

A comment may be in order: although we have discussed only the 3-point 
function in this section, it is clear that the same method as above can be applied 
to calculate arbitrary n-point Green functions of current operators. Namely, we can 
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simply close the open fermion legs of the unamputated vertex function G^^'^^^~^ 
by inserting an n-th vertex factor 7^^". The reader can apply this procedure, for 
instance, to T^'^^ given in Fig. 9 in the ladder approximation and can see which 
types of the diagrams give gauge invariant 4-point Green function for photons 
corresponding to 77 — > 77 scattering. This example shows that even in the ladder 
approximation there generally appear corrections intrinsic to the Green functions 
which can be attributed neither to the propagator nor to the gauge boson vertices. 



5. Low Energy Theorem and Pseudoscalar Bound State 

Thanks to these properties of our propagator and vertices, the decay amplitude 
for 7r° 27, for instance, satisfies the low-energy theorem. This holds, of course, 
provided that we use BS amphtude for the boundstate 7r° calculated by the 'same' 
approximation as that used for propagator and vertices. To see this is the purpose 
of this section. 

We consider the 3-point function T^^^ip^k) for on-shell photon case, — 
— 0, and recall the following fact which holds in the full theory and in our 
approximations at any order also. The vector gauge invariance {Pfj,T°'^'^{p, k) — 
ki,T'^'"'{p,k) = 0) and the bose symmetry T'^'"'{p,k) = T'^''t'{k,p) together with 
odd parity property, determine the most general form of T°^i^'^ as* 

T«^'^(p, k) = e^'^P^Ppkaq^'Fiiq^) + {e'^'P^k'' - e'"'P^pP) PpkaF2{q^) . (5.1) 
★ Another possible form 

can actually be expressed in terms of the above two amplitudes: Indeed, from the identity 

gliJ^'^P'^qa] = 0, we can prove 

A"^"'{p, k) = e'^'P'^Ppkaq" - {e'^^'P^k" - e^^P'^p'') Ppk„ . 



23 



Therefore generally we have 



gaT^^'^lp, k) = e^^-Pypk,q^F^{q^) . (5.2) 

The explicit factor in the RHS implies that only the pole term in -Fi(g^) 
can contribute to this amplitude qaT'^^^ in the limit of on-shell pion, q^ 0. 
(This is the reason why there exits a low-energy theorem to the tt^ 27 decay 
amplitude.) It is important to note that the vector gauge invariance is essential 
here to this conclusion; indeed, otherwise, the other invariant amplitude term of 
the form €'^^^'^P{pp — kp)Fs{q'^), for instance, could appear in Eq.(5.1) and contribute 
to this amplitude at = without having massless pole. 

We now also recall the definition of the BS amplitude: the BS amphtude x as 
well as its conjugate x for the pion boundstate |7r(q)) normalized invariantly as 
(7r(q)|7r(q')) = (27r)32|q|53(q - q'), are defined by 



(0| Ti;{x)i;{y) |7r(q)) = e'^^^ J ^ e-'P' x{p; O) , 
(7r(q)| TV^(y)V^(a:) |0) = e^"^-^ j ^ e*^^ x(p; ?) , 



(5.3) 



where X = {x + y)/2 and r = x — y. It is easy to see that CPT invariance leads to 
the relation x{p] q) = x{P] because of which we often denote x as x below for 
notational simplicity. This boundstate generally appears as a massless pole term 
in the 4-point amplitude G4 = (0 1 Tipipipip \ 0) in the form: 



G4 = X -9 X + regular term at g'^ = . (5.4) 



If we close the first -0 and -0 legs of G4 by inserting ^7") it becomes the axial- vector 
vertex with fermion propagators attached. Therefore the pion pole appears in 
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the axial- vector vertex in the form (see Fig. 15): 

= i /tt^^^X + (regular term) , (5.5) 

where x is the amputated BS amplitude, x = S^^x^f^^ decay constant 

/tt is given by 



i/.g" = (0|VS(0)y7"75V'(0)|7r(q)) 




Fig. 15. BS amplitude of the pion, x, in the axial-vector vertex, Fg . 

These are general arguments. If we use an approximation to determine G4 and 
F5, the BS amplitude x is also determined accordingly, via Eq.(5.4) or (5.5): for 
instance, if we use the ladder approximation (/C2PI = ^2Pl)i then the 4-point 
function G4 is the ladder L given in Fig. 6 and therefore the BS amplitude x should 
be determined by the homogeneous ladder BS equation as shown in Fig. 16. The 
decay constant /tt is also the one calculated by Eq.(5.6) using this BS amplitude. 




Fig. 16. Graphical representation of homogeneous ladder BS equa- 
tion for X. 
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The massless pole in T"^^ required in Eq.(5.2) is provided by the pion pole 
in the axial- vector vertex as given in (5.5). In the ladder approximation, the 
3-point function T'^^'^ is given in the form (4.6) so that we find 

lim qaT'^^^p, k) = ie^TrTr {xS^r S^T^" S^) + (cross term) , (5.7) 

which is shown in Fig. 17. 

P 




k 

Fig. 17. Graphical representation of the function T^'^(p, /c). 

Combining with Eq.(4.14), this proves that 

r'''^(p, k) = ie^Tr (x-SFF'^^pr^-Sp) + (cross term) = -^[Anomaly] . (5.8) 

Namely, when the BS amplitude x is calculated in the ladder approximation, then 
the 7r° — > 27 amplitude consistent with the low energy theorem can be obtained 
by the LHS of Eq.(5.8) with vector vertices and fermion propagator which are also 
calculated by the same order of approximation. It will also be easy to convince 
oneself that in the next order of approximation taking /C2PI = ^2PI ^" ^2PI^' 
correct amplitude satisfying the low-energy theorem can be obtained by calculating 
the diagrams given in the RHS of Fig. 13 with Fg vertex replaced by the BS 
amplitude x- 
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6. Saturation by zeroth order term 



In the ladder approximation, there is a gauge in which the solution to the SD 
equation take the form 

i5Fi(p)=y-E(p2) (6.1) 

without wave function renormalization. [This is the case in Landau gauge for fixed 
coupling case*, and there is a non-local gauge realizing Eq.(6.1) even for running 
coupling case [3,10].] In such a case, there appears an interesting phenomenon 
that the saturation of the low-energy theorem by the LHS of Eq.(5.8) is already 
realized by the zeroth order terms in the eocternal momenta of the BS amplitude x 
and the vector vertices and . This is actually surprising because the linear 
terms in the external momenta individually give non-vanishing contributions to 
the amplitude. The net effect of those linear terms vanishes. Let us explain this 
fact. 

First let us see what is the zeroth order term in the BS amplitude x- The WT 
identity Eq.(3.9) for the axial- vector vertex reads 

,.r?(p.|.p-|)^..-(p.|)|-..n4s.'(p-|) . (.2) 

Consider — limit on both sides: For the LHS, only the pseudoscalar pole 
term given in Eq.(5.5) survives in this limit and yields 



lim LHS = -Ax(p;g = 0) . (6.3) 

g->0 



When S-p takes the form Eq.(6.1), the RHS, on the other hand, gives 



lim RHS = -2E(p2)^75 . (6.4) 

q^O 2 



★ If we use the regularized gluon propagator, (2.4), there is a shght deviation from Landau 
gauge in the region ^ [3] . 
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So we find 

g = 0) = -^E(p2)^^5 = ;^o(p) . (6.5) 

In ^ 

This is an exact result giving the zeroth order term in q as far as S'p takes the 
form Eq.(6.1). [Pagels-Stokar's proposal [11] is essentially to regard this as an 
approximate BS amplitude up to order g^: 

As a matter of fact, there appear three terms which generally take the form 

X(P, 9) =7-S(p^)^75 + Xlinear + C'(g^) , 

. . 1 X (6-7) 

Xlineax = ( " ?¥+ Qiv'M^ T{v\W- ) 75 ■ 

So the Pagels-Stokar approximation corresponds to discarding all those am- 
plitudes.] 

Next let us see the zeroth order form in the vector vertex F^. Again from the 
vector WT identity (3.9) and the form (6.1) of S^-, we find 



v^^^ + 1' ^ - f ) = ^"^F' + 1) Q - Qis^^ - f ) 

with 

2fc+l 



(6.^ 



k=0 

for p2 — (on-shell photon) . This tells us that 



r'^ + f , ^ - 1) = - 2e^f{e^p- 1) + rf,] g , (e.io) 

where F^j. stands for transversal part satisfying Pj^F^j. = by itself. Since 
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generally takes the form 



+ (p> _ p^^p . [Ci{f,p- 1) +/C2(£2,p • I) + /C3(£2,p • l)\ (6-11) 

and is at least linear in the external momentum p, the zeroth order term in p of 
r'^ is given by 

f£ + -, £ - -) = 7^ - = T^{£) . (6.12) 

V 2 2/ p=o 

Note that here stands for the average momentum of those of incoming and 
outgoing fermions. 

Let us now show that the zeroth order terms xo a-nd Fq of the BS amplitude and 
the vector vertices already saturate the low-energy theorem; namely, the equality 

1 



T''^{p,k) =%^\p.k) = f [Anomaly], 



(6.13) 



T^^^ip, k) = ie^Tr {xoSf^oSf^qSf) + (cross term) (6.14) 

holds for the on-shell photons and pion (p^ — k'^ — = Q). A key observation for 
this equality is that 



'^3 a 
5a ^7 75 



■ S(£ + pj) ^75 + 75^ - S(£ - k)) 



+ (S(£ + p) + E(£-A;))|75 



^Sp\i + p)^^5 + jl5tS^\i - k) + xo [i + 2 



T3 



p — k 



(6.15) 



Eliminating xo using this identity, we find 



+ 



Tr (i^^^SFr^SF^) + Tr (5F^^5F^|75^S) 

IV (^Fi|75r^^Fr^) + Tv (^Fr^i|75^Frs) 

QaTr ^^7"75S'FroS'Fr(^5'F^ + (cross term) 



(6.16) 
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since the last 0{q^) term in Eq.(6.15) cannot contribute to the on-shell amphtude 
after the loop integration. Although we have suppressed the momentum arguments 
here, it is important to note that the factor ^75 indicates also the place at which 
the momentum q flows in, so that it is not anticommutative with the vector vertex 
(or Fq) which has nontrivial dependence of the leg momenta. The LHS is well- 
defined since the loop integration is convergent because of the presence of the BS 
amplitude xo which damps sufficiently rapidly. But the terms in the RHS are not 
well-defined separately. In order to make each term well-defined, we add to the 
RHS the identity for the PVG propagator 5f, 

= -^e2(^5Ffey7"755F - i^i^S^ - Spt^^^ - 2MSf^^5Sf)QYSfQ7^ 
(following from Eq.(4.12) ) and its /i u exchanged one. Then Eq.(6.16) reads: 



le 



Tr (t^^^SF^SF^) + Tr (SFT^SFt^^sn) 
V 2 ^ /PVG V ^ 2 /PVG 

Tr (SFt^^sT^SFn) + Tr (SFT^t^^sSFn) 

\ 2 ^ /PVG \ "2 /PVG 

g«Tr (^7«755Fr^SFr(;5F) + (cross term)] 

2MTr (y75'5FQ7z.'5FQ7/x'5F) + (cross term)] , 



-|- ie lim 

M— »-oo 

(6.17) 

where each trace term with suffix PVG denotes the regularized one obtained by 
subtracting the same form term with replacements Sf — > Sf, Tq^'^^ — > Q^^^^^ 
made. For instance, the first term reads 

Tr (t^^^SFV^^SFvi)^^^ = Tr (z|755fF^5fF[^ - z^^^SfQ^SfQ^) ■ 

(6.18) 

In this regularized form we can make the shift of loop integration variables freely 
{i.e., producing no surface term) as explained before. Then the first four trace 
terms in Eq.(6.17) cancel among themselves and vanish: Indeed, in order for such 
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'biangle' diagrams with 75 inserted to give non-vanishing contribution, at least 
four 7 matrices have to be picked up, so that only the (momentum-independent) 
elementary vertex parts (^7^ and (57^^ in Fq and Tq can contribute. But then, the 
factors 1^15 anti-commute with Fq or Fq and those diagrams cancel with each 
other. This is the case since the shift of integration variables is allowed by the 
above 'regularization'. But, this 'regularization' introduced an extra contribution 
given by the last term in Eq.(6.17), which is precisely the anomaly required by the 
low-energy theorem. So what now remains to be proved is that the second term in 
Eq.(6.17) vanishes: 

T^^'^ip, k) = i^e^Tr (^^^""^^Sf^^SfT^Sf) + (cross term) . ^^'"^^^ 

This can be proven if we can show that the amphtude Tq^'^ satisfies (vector) gauge 
invariance, p^T^^^ — kuT^^^ — 0. Indeed, then, together with the Bose symmetry, 
Tq^^ has to have the general form Eq.(5.1) and hence QaT^^^ can be nonzero at 
= only when a massless pole term exists in the axial- vector channel. 
But the diagram of Tq^^^ is clearly regular at = since the fermion here carries 
non-zero mass function S, and hence QuTq'^'^ = on the mass shell = 0. 

To show the gauge- invariance of Tq^'^, we note that we can replace the zeroth 
order vertex function Fq and Fq in Eq.(6.19) by the full non-transverse parts 7^^ — 
2^^'f{f,p ■ £) and -f" - k-l), respectively, defined in Eq.(6.8), since the 

differences are of e)(p2) and 0{P) and vanish on mass shell {q^ — p"^ — k^ — 0) 
after loop integration. Then, with these replacements performed, we can use the 
vector WT identity Eq.(6.8) to obtain 



pX^'ip, k) = i^e^ [Tr (^7"75^Fr^5FiQ) - TV {^-^iQi'^i^SfTISf 

Tr (|7"75^F^^^Q^F) - Tr {^-^^'^^^SfiQTISf) 



(6.20) 

Here Q indicates also the position from which the momentum p flows out. The 
first two terms clearly cancel with each other. The second two terms also cancel 
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since, again, only the (momentum-independent) elementary vertex part Q^'^ in Fq 
can contribute because of the presence of 75 in the trace and hence the factor Q, 
indicating also the momentum-flow position, becomes commutative with Fq. This 
finish the proof. 

It is interesting to see explicitly how the diagram (6.14) with lowest order terms 
Xo, Tq and Fq inserted actually reproduce the value of the low-energy theorem: 



i(2^ 



tr 



+ 



P 



k 



1 



S(£ - k) 



i n 



+ (p ^ /c, /i. z/) 




2J S(£ + p) 



xE(x) (E(x) - 2xE'(x)) 



47r^ fn J (E2(x)+x) 



(6.21) 

where x — —i'^ denotes the Euclidean momentum. Interestingly enough, the value 
of the integral in the last line is ^ irrespectively of the detailed functional form of 
(nonzero) E(x). This can easily be seen if one can make a change of the integration 
variable x into dimensionless one y — x/Tp'ix). Thus the lowest order terms 
correctly reproduce the value of the low-energy theorem. 
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APPENDIX A. Dimensional Regularization 



If we use the dimensional regularization instead of PVG regularization, the 
identity used in the algebraic manipulation deriving the WT identity Eq.(4.5) is 
changed to 




where i. is the (now n-dimensional) loop momentum and / is its extra dimensional 
part {^ — + I) . The first 3 terms in the RHS is the usual one and the usual 
proof of the WT identity applies to those terms: then, they vanish since now the 
integral is regularized and the shift of the integration variable is justified. So we 
have to evaluate only the contribution from the last term coming from the extra 
dimension part I: 

qaT''^"'{p, k) = 2ie^laTr (r^^Fr^^pr^^p) + (cross term) . (A.2) 

Since I oc n — A, only the divergent parts can contribute in the RHS. Since the 
propagator S-p and the vertices F'^, and Fg asymptotically behave as free one 
as explained in (4.9), we can evaluate as follows: 

(RHS) = hm 2ze^Tr (/75'5fs=oQ7z.'5fe=oQ7/.'5fe=o) • (A.3) 

n— »-4 

But this is again the same diagram which leads to the well-known anomaly. 
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